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Введение. Увеличение скоростей движения железнодорожного транспорта и повышение нагрузок на оси колес- 
ных пар обуславливают необходимость модернизации существующего парка. Научные исследования в области 
динамики подвижного состава направлены на учёт колебательных процессов, возникающих при движении же- 
лезнодорожных транспортных средств в традиционном конструктивном исполнении. Присоединение дополни- 
тельных элементов рассматривалось на уровне сцепки двух вагонов и присоединении третьей тележки в центре 
тяжести железнодорожной платформы. Построению математических моделей, позволяющих оценить динамиче- 
ские состояния таких конструктивных решений, в научной литературе не уделено достаточно внимания. Цель 
данного исследования — создать метод оценки динамических состояний вагона. Рассматривается ситуация, ко- 
гда в его структуру вводится дополнительная совокупность масс-инерционных и упругих элементов, причем от 
корректировки их параметров зависит общее динамическое состояние транспортного средства. 

Материалы и методы. Базовым инструментом проведения исследований является структурное математическое 
моделирование, в основе которого лежит подход, когда исходная расчетная схема представляет собой механиче- 
скую колебательную систему в виде твердого тела на упругих опорах с дополнительной введёнными в её струк- 
туру типовыми элементами. Динамическим аналогом используемой расчетной схемы является структурная схема 
системы автоматического управления, применение которой позволяет детализировать связи между типовыми 
упругими и масс-инерционными элементами. 

Результаты исследования. Предложен метод оценки динамических состояний железнодорожных транспортных 
средств, основанный на построении математических моделей, с учетом введения дополнительной структуры 
масс-инерционных и упругих элементов. Исследовано влияние дополнительных параметров на динамическое 
состояние транспортного средства. Получены аналитические соотношения, позволяющие при изменении соот- 
ветствующих параметров технического объекта снизить динамические нагрузки на основные конструктивные 
упругие элементы. Приведена передаточная функция межпарциальных связей, позволяющая контролировать вза- 
имодействие между координатами движения транспортного средства при действии двух кинематических возму- 
щений синфазного типа. 

Обсуждение и заключение. Сформированная математическая модель позволяет оценить динамическое состоя- 
ние железнодорожного транспортного средства в условиях действия кинематических возмущений. Результаты 
исследований могут быть применены при модернизации существующих и создании новых транспортных средств 
с улучшенной динамикой. 
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моделирование, структурная схема 
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Введение. Расширение технико-экономических связей межрегионального уровня, обеспечение роста произ- 
водственно-технического потенциала страны и поддержание разрастающейся системы международных торгово- 
коммерческих отношений во многом зависит и опирается на железнодорожный транспорт [1]. Объёмы перевозок, 
выполняемых железнодорожным транспортом, постоянно возрастают, что обусловливает необходимость учета 
нежелательных эффектов от повышения динамических нагрузок [2], напрямую влияющих на надежность эксплу- 
атации как подвижного состава, так и верхнего строения пути. Несмотря на наличие таких негативных факторов, 
необходимо выполнять плановые показатели, к которым можно отнести повышение участковой скорости, со- 
блюдение норм веса и длины составов, увеличение осевой нагрузки до 30 тонн и более. Это отражает реальные 
запросы развития российской экономики и стимулирует создание новых более мощных локомотивов, обновление 
парка подвижных средств, модернизацию путевого хозяйства [3]. Вместе с тем нельзя не рассматривать и воз- 
можности проявления негативных последствий интенсификации перевозочных процессов. Одним из важных во- 
просов является возрастание или увеличение темпов износа верхнего строения пути с соответствующими выте- 
кающими сложностями [4]. 
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В текущий момент времени повышенное внимание уделяется развитию методологии оценки динамического со- 
стояния подвижного состава, взаимодействию технических средств с рельсовыми путями, экономии электроэнер- 
гии, повышению надёжности и безопасности перевозочных процессов. Методология математического моделирова- 
ния описана, например, в [5]. Вместе с тем, существуют и другие возможности для поиска рациональных реше- 
ний [6]. Важно уделить внимание вопросам модернизации существующего парка грузовых вагонов, эксплуатация 
которых при повышенных нагрузках уже не является эффективной [7]. Одним из подходов, который мог бы быть 
принят к разработке, является концепция установки для грузовых 4-осных вагонов дополнительной двухосной те- 
лежки [8]. В этом случае можно ожидать более равномерного распределения нагрузки на верхнее строение пути, а 
также возможность увеличения веса перевозимых грузов при сохранении регламента на значения осевой нагрузки 
в пределах 22 т [9]. Неравномерность параметров контакта колесо-рельс инициирует колебательные движения ва- 
гона, что формирует в свою очередь колебательные движения корпуса вагона. Процесс колебаний формируется 
также условиями взаимодействия вагонов внутри состава [10]. В этом случае возникают динамические реакции 
связей, накладывающиеся на статические составляющие полной реакции, что может существенно увеличить уро- 
вень динамических взаимодействий в контакте колесо-рельс [11, 12]. Однако возможностям структурного матема- 
тического моделирования при оценке динамических состояний железнодорожных транспортных средств при вве- 
дении дополнительных связей пока не уделялось достаточно внимания. Поэтому целью данного исследования яви- 
лось формирование метода оценки динамических состояний железнодорожного транспортного средства при введе- 
нии в его структуру дополнительной совокупности масс-инерционных и упругих элементов, корректировка пара- 
метров которых будет влиять на общее динамическое состояние транспортного средства. 

Материалы и методы. Методологической основой исследований выбрана структурная теория виброзащит- 
ных систем, позволяющая в линейной постановке с учетом сосредоточенных параметров и малых колебаний от- 
носительно положения статического равновесия или установившегося процесса, достаточно точно оценить ди- 
намические свойства железнодорожного транспортного средства. Расчётной схемой является механическая ко- 
лебательная система с динамическим эквивалентом в виде структурной схемы системы автоматического управ- 
ления. Это позволяет детализировать связи между элементами транспортного средства, а также использовать 
методы, характерные для теории автоматического управления (передаточные функции, преобразования струк- 
турных схем, свертки и упрощения, частотные характеристики) [13]. 

Стандартные конструктивные схемы грузовых железнодорожных транспортных средств описываются извест- 
ными расчётными схемами, а их динамические особенности можно оценить с использованием линейных расчёт- 
ных схем. Рассматривается железнодорожное транспортное средство в виде грузового четырехосного вагона, 
предназначенного для перевозки тяжёлых грузов, таких как каменный уголь, руда, песок, металлопрокат, мало- 
габаритные металлические конструкции и др. (рис. 1). 
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Рис. 1. Принципиальная схема грузового четырехосного вагона 


Структура представленного железнодорожного транспортного средства содержит кузов, обладающий мас- 
сой М, моментом инерции Ји опирающийся на две двухосные тележки, условно представленные как совокупно- 
сти масс-инерционных и упругих элементов. Динамические особенности рассматриваемой системы показывают 
наличие избыточного влияния на конструктивные элементы тележек железнодорожного транспортного средства. 
В этой связи для улучшения его динамических СВОЙСТВ В структуру подвески дополнительно вводится третья 
двухосная тележка, которая размещается в центре вагона (рис. 2). 
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Рис. 2. Принципиальная схема грузового шестиосного вагона 


При собственном весе вагона около 23 тонн и осевой нагрузке в 22 тонны, это может обеспечить увеличение 
веса перевозимых грузов ориентировочно на 20 тонн, то есть существенно повысить эффективность перевозоч- 
ных процессов, не создавая излишних динамических нагрузок на верхнее строение пути (ВСП). Практическая 
реализация предлагаемого подхода требует модернизации конструкции типового вагона, которая заключается в 
создании узла крепления дополнительной тележки для уменьшения нагрузки на ось грузового вагона за счёт 
введения дополнительной двухосной тележки. Модернизация узла крепления осуществляется в тех же конструк- 
тивно-технических формах, что и крепления с использованием шкворней в двух «штатных» двухосных тележках, 
то есть через установку надшкворневой и подшкворневой балок с соответствующим шкворневым узлом стан- 
дартного типа. 

Предлагаемая методика повышения эффективности использования четырехосных грузовых вагонов в усло- 
виях дополнительных нагружений ориентирована на решение проблемы повышения нагрузки на ось колёсной 
пары до 30 тонн и увеличении скорости движения поездов. Это достигается модернизацией типового грузового 
вагона путём установки дополнительной двухосной тележки соответствующим устройством, обеспечивающим 
условия её динамического взаимодействия с конструкцией рамы грузового вагона. 

Установка дополнительной двухосной тележки за счёт перераспределения нагрузки между общим набором 
колёсных пар обеспечивает возможности перевозки грузов больших по весу при одновременном снижении 
нагрузки на ось. Это обеспечивает более рациональные условия эксплуатации рельсового пути и верхнего стро- 
ения железнодорожного полотна при сохранении приемлемой длины состава. 

Предполагается, что шкворневый узел будет иметь упругую резиновую прокладку, обеспечивающую амор- 
тизацию взаимодействия тележки и рамы корпуса вагона. При этом не предполагается модернизация дисков ко- 
лесных пар. Изменяются лишь форма проводки элементов «штатной» пневматической системы торможения со- 
става и конфигурация проводки пневмопроводов. 

Результаты исследования. Расчётная схема рассматриваемого железнодорожного транспортного средства в 
первом приближении может быть представлена в виде механической колебательной системы, состоящей из твёр- 
дого тела, обладающего массой т, и моментом инерции 7, опирающемся на три упругих элемента с жесткостями 1, 


Ко, о. Кинематические воздействия представлены синфазными гармоническими функциями одной частоты (рис. 3). 


У М,Ј 


21 


Рис. 3. Расчётная схема железнодорожного транспортного средства 
при кинематическом возмущении (21, 22, 20) 
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Центр масс системы — точка О находится на расстояниях / и р от концов твёрдого тела (точек 41, В). Эле- 
мент с жесткостью Ао закреплён в точках Р и 1; расстояние ОР. обозначено как /о. Движение системы рассмат- 
ривается в координатах ут, уз и уо, ф, связанных с неподвижным базисом. В расчётах используются соотношения: 

Уо = ауу +Бу›, ф = с(у› ул), у= уо 1ф, уз = уо +, 
1 1 1 (1) 
Э 1 
Ур = уо +1ф, а = „Б = ‚с = . 
1+1, 1+1, 1+1, 


Для вывода дифференциальных уравнений движения используется формализм Лагранжа [14], что требует по- 
строения выражений для кинетической и потенциальной энергий. В данном случае имеем: 


2 


г = т(ау кру) +? (у-у), (2) 


1 1 1 
П= № (у а)? + 5 (Ур 3) 02 2): (3) 


Учтём, что потенциальная энергия (3) может быть записана также в виде: 


1 1 1 
П Б (у -21) Е (> -2,) Ы С? + руз +1с(у2 -у)-2] Е 


1 2 1 Б 2 (0 
А (2 –21) "25 (уг —22) Ко [в + у –2, | А 
где а = а – [0с, Б = Б + ос. 
Система уравнений движения в координатах у, у во временной области примет вид: 
я (та? + Је? )+ У (А + ка?) У 2 –таһ) + уз Коа = 121 + Коа120, (5) 
уз (ть? + с? + у» (к +?) у (с? —тав) + уа = 2 + Кова. (6) 


После применения интегральных преобразований Лапласа при нулевых начальных условиях [15] система 
уравнений (5), (6) может быть представлена в операторной форме: 


У (та? + Је? )р? + +? |- У» (с> — тар) р? - ка | = 2 + 00120, (7) 
ӯ, (С +?) р? +8, +60 Е (с — таБ) р? — Киа | = 2, + БЫ, (8) 


где р = јо — комплексная переменная ( / = \/—1 ), значок <> над переменной означает её изображение по Лапласу [5]. 

На основе (7), (8) построена, в соответствии с положениями метода структурного математического модели- 
рования [5], структурная математическая модель в виде структурной схемы эквивалентной в динамическом от- 
ношении системе автоматического управления (рис. 4). 


(Је? – Ма) р* -Кав, 


1 1 \ 


2_ 2_ "9 
^ (Ма? + Јс?)р? +в + Коа? ао ае (МБ? + Је?) р? + + Корр У 
а я Е 
а «2 РР 88 
К а /% 2 | 20 
061 Кор! = 


Рис. 4. Структурная математическая модель железнодорожного транспортного средства 


Обсуждение и заключение. Особенности системы заключаются в том, что она имеет два парциальных блока, 
каждый из которых определяет соответствующие парциальные частоты: 


2 
2 _ А + Коа! 


Е, 9 

} та? + Јс? 9) 
к + кор 

Э.а. 2 0-1 

7? > тр? + 2’ (10) 


которые, в свою очередь, определяют границы расположения частот собственных колебаний системы в целом: 
2 2 2,2 2 
Особ < И < 05 < 05, (11) 
ГДЕ б91соб, @2соб — частоты собственных колебаний системы, при работе на которых возможно возникновение ре- 
зонансных режимов. 
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К числу особенностей системы относится наличие нескольких одновременно действующих внешних возму- 
щений. Полагая в целях упрощения, чтот 2, = 2, = 2, (вполне допустимо на стадиях предварительных динами- 


ческих оценок) примем, что на вход первого и второго парциальных блоков действуют силовые факторы: 
О = 2(8 +Коа, ), (12) 
О, =2 (К +). (13) 
Используя структурную схему нарис. 4, запишем выражения для передаточных функций, полагая, что между 
внешними факторами возмущения имеется связь, формируемая соотношением: 


О, = 0, 0, = оО, (14) 
Риа (к + а) (тб? +?) р? +, + + (6 +в, (2° та) р? сһаһ). Р 
2 А(р) 
я, ()- У: _ о (5 +5, | (та? +? р? + + а? |+ (6 + а (ое? -таһ)р? ња) о 
2 А(р) 
где 
А(р) = (та> +?) р? + +Коа? | (ть? +?) р? + +? |-[ (2? — таЬ) р? аһ |, (17) 


является частотным характеристическим уравнением системы. 
Числителем передаточных функций в выражениях (15), (16) определяют режимы динамического гашения ко- 
лебаний, что может быть детализировано из уравнений, получаемых при «обнулении» числителей (15), (16): 
2 
р (А + оа) (К + об? )- а. (6 + № Кое 8) 
09 |дин т Й 
(А + а; ) (ть? + Је? ) +а( + Кор)(Је? – таь) 


К (о + ) (в + а? )- о (А + оа) оар, у 
(А +6, (та? + Ле?) + (№ + аи (с? –таБ) 


Из выражений (18), (19) следует, что в системе с двумя степенями свободы при наличии двух связанных 
между собой возмущающих факторов возможно возникновение режимов динамического гашения колебаний на 
двух частотах, параметры которых зависят от значений коэффициента связности а. Этот коэффициент может 
принимать отрицательные, положительные и нулевые значения. 

Из анализа структурной математической модели (рис. 4) следует также и возможность возникновения особого 
динамического режима на частоте: 

Е Ко 

"Л? тар’ 
когда парциальные блоки получают возможность разъединения. В этом случае система (рис. 3) распадается на 
два автономных блока, которые не будут создавать ситуаций взаимодействия парциальных структур. 

Реализация такого режима может привести к существенной разнице в значениях отклонений в точках Ат и В! 
и «разбросу» значений динамических реакций в точках А1, Ву и О1. Определение динамических реакций в точках 
А, Вир может быть реализовано по методике, изложенной в работе [14], в которой динамическая реакция опре- 
деляется как произведение динамического смещения (например, точках 41, Ву и Р)) на значение приведённой 


(20) 


динамической жёсткости. 
Для координаты у, динамическое смещение определяется из выражения для передаточной функции (15), а 


для координаты у,, из выражения (16). Для определения приведённой динамической жёсткости используется 
частотное характеристическое уравнение (17): 


Є –тађ)р? -каһ | 


Коа?а! 
Юн (р)= + – . (21) 
я ( (ть? +?) р? + +2 (ть? +?) р? + +2 
Аналогично может быть найдено значение приведённой динамической жёсткости по координате у, 
2 
2 2 
Е Б?Ь (Је — та) р — Куа 
про (р) = К в о . | | (22) 


(та? +?) р? + + оа? (та? +?) р? + + оа? | 
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Для определения конкретных значений К. пр ( р) С ( р) необходимо найти значение частоты 6, которая 
определяется из условия, что у, / у, =1. В более общем случае предполагается, что: 
22 1, (23) 
У 
где ү — коэффициент связности движения элементов по координатам у, и у;. Таким образом этот коэффициент 
(для случая у = 1) может быть записан в виде: 
у, @ (5, + Кор; (та? + Л? )р? + + Коа? | + (А, + Коа (ое? – тар) р?- Коа | 
У (А + Коа} („> + Је? )р? + +02 |+ о (5 + Кору )] (<? —таБ) р? ар | 


(24) 


Принимая конкретные значения у, можно найти частоты колебаний для движения рассматриваемого объекта 
в координатах у, и у. Например, при ү = 1 частота поступательных вертикальных колебаний твёрдого тела опре- 
делится выражением: 
у а (+в ) (к + оа? ) (& + оа) Ков + (Ки + Коа ) коа, — 
я а (№5 + )(та? +? )+(, + Коа ) (<? та) 
(А + оа) (6 + об?) (К + Ко ) (Ки + оа?) +9 (8 + Б) а 


(5 + а ) (тБ? +? )-в(®, + ) (<? — таБ) 


(6) 


(25) 


Частота © синфазных гармонических возмущений 21(0, 22(1) и 20(1) обеспечивает движение твёрдого тела при 


Фф = 0, то есть у = У, = ур. При других значениях у (у = 1) определяются у, и у,, а на их основе из геометриче- 
ских соображений определяется значение ур. Для твёрдого тела при известных у! и у, легко может быть 


найдено положение центра вращения (или колебаний) [12]. 


Динамические реакции связей Ад, Кр и Кр, могут быть в первом приближении найдены по формулам: 
Ка = оу, Ки = КУ, Кр = Коур. 
В общем случае динамические реакции определяются с использованием динамических смещений у, У, ур, 


определяемых выражениями (15), (16). Что касается динамического смещения в точке О, то используется выра- 


жение ур =ау+Ру›, параметры которого определяются вышеприведёнными значениями. В выраже- 


ния (15), (16) могут быть введены данные о связности силовых факторов воздействия (параметр о). Для получе- 
ния конкретных данных о значениях динамических реакций вводятся параметры частоты, при которой реализу- 


ется соотношение амплитуд колебаний у, и у, через коэффициент связности амплитуд колебаний ү. 


Полная реакция связей в точках А, В и Р определяется суммой двух компонент: статической и динамической. 
Статическая компонента может быть найдена из выражения для передаточных функций динамических реакций 
при р = 0 и задании параметров статической нагрузки (вес вагона и перевозимого груза). При интенсивном раз- 


витии колебательных процессов, когда возрастают колебания по координатам Уь Уз А Ур, полная реакция может 


изменяться в значительных пределах и отличаться от статической реакции связей. При наличии динамической 
компоненты полная реакция может принимать различные значения, в частности, нулевые или отрицательные. 

Сформированная в рамках предложенного метода математическая модель, обозначенная выражением (25), 
позволяет оценивать динамическое состояние железнодорожных транспортных средств при введении дополни- 
тельных связей в их структуру для формирования комплекса рекомендаций по получению устойчивых режимов 
их эксплуатации. Изучение особенностей системы при помощи подходов, характерных для структурного мате- 
матического моделирования, позволяет детализировано рассмотреть связи между элементами. Применительно к 
рассматриваемому техническому объекту это даёт возможность корректировки динамического состояния техни- 
ческого объекта на основе варьирования параметров совокупности дополнительно введённых элементов и поз- 
воляет уменьшить нагрузку на основные части подвески, а также установить наличие в системе собственных 
частот и частот динамического гашения колебаний. 
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В дальнейшем планируется проведение исследований при введении в структуру транспортного средства 
демпферов и устройств преобразования движения для оценки возможностей структурного математического мо- 
делирования. Также интересным направлением является оценка возможностей изменения динамических реакций 
в зависимости от внешнего воздействия, что позволит оценить усилия, прилагаемые к различным элементами 
подвески транспортного средства. 
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